Tunable anisotropic nonlinearity in superconductors with asymmetric antidot array 
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The influence of the spatial asymmetry of the pinning potential on the spectral composition of 
the voltage, induced in perforated superconducting Al bridges by the injection of a sinusoidal bias 
current, was investigated. The loss of the mirror symmetry of the pinning potential leads to the 
appearance of even Fourier components in the induced voltage in the vicinity of the superconducting 
phase transition line on the H — T diagram [H is the external magnetic fleld, T is temperature). 
Artificially-introduced asymmetry for vortex motion makes it possible to create low-resistive mate- 
rials, in which nonlinearity depends on the direction of injected electrical currents. 

PACS numbers: 



The controlled manipulation of magnetic flux quanta 
(or vortices) in superconducting films under the action of 
non-equilibrium fluctuations or external excitations with 
zero average (vortex ratchet) has recently attracted a lot 
of interests (see, e.g. , Refs. % i, i, i, % i, 7, 1, J)., E [u 
121 . Il3l . [ij, ll5| and references therein) . A preferable vor- 



tex motion in a certain ('easy') direction can be exper- 
imentally realized in various one-dimensional (ID) and 
two-dimensional (2D) systems of reduced symmetry: in 
films modulated by triangular dots 0] , in perforated films 
with regular arrays of asymmetric antidots lol . ITol . 11 1 , 
in mesoscopic singly-connected superconductors [l2| , and 
in superconductors placed in a spatially-modulated mag- 
netic field [13, 13, Il5l |. All these systems can be poten- 
tially used for designing flux pumps for removal of un- 
wanted trapped flux in superconducting devices [l[, and 
dc rectifiers (diodes). 

In this Letter we are aiming at the generalization 
of the concept of unidirectional vortex dynamics in a 
ratchet potential for the case of higher harmonics gener- 
ation. Indeed, passing an alternating current of certain 
frequency / through any medium should generally in- 
duce voltage oscillations at multiple frequencies nf (n 
is integer). However, the appearance of even harmon- 
ics (regardless on the frequency range) is known to be 
forbidden for a material with spatial inversion or mirror 
plane symmetry [l^ . One can expect that an artificially- 
introduced asymmetry for vortex motion, occurring in 
all ratchet systems, will lead to a nontrivial nonlinear 
response. In the presence of a ratchet potential the 
static I — V dependence of the perforated superconduct- 
ing bridge will lose symmetry, V{I) ^ —V{—I). Since 
the voltage drop in a sample should be zero at / = 0, the 
expansion of V vs. / in power series has the following 
form: V — ail + a2l^ + a^P + where are con- 
stants depending on the external magnetic field H and 
temperature T. Substituting I{t) = Jo sin(27r/<), one 



can see that the time-dependent voltage V{t) oscillates 
at multiple frequencies, including zero frequency. Pro- 
vided V{I) = —V{—I), both the diode effect and the 
even Fourier components should disappear. It is worth 
noting that breaking the sample symmetry in 2D super- 
conducting films can affect only a certain direction, while 
keeping the symmetry in the perpendicular direction un- 
altered. It gives us a possibility to fabricate an artificial 
material with anisotropic nonlinearity, depending on the 
orientation of the injected electrical currents. 

Two perforated Al microbridges of the width 600 ^m 
and thickness d = 50 nm, which have the critical temper- 
atures Tc close to 1.3 K, were prepared by dc-sputtering 
on a Si/Si02 substrate. The pinning landscape, pre- 
pared with electron beam lithography, has the period 
a = 3000 nm and it can be represented as a superpo- 
sition of two interpenetrated square sublattices of big 




FIG. 1: (Color online) (a) A schematic view of the perfo- 
rated superconducting Al bridge; (b) and (c) Atomic force 
microscopy images of samples A and S. White dashed lines 
depict the unit cells (3x3 /iia^). 
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FIG. 2: (Color online) An example of the non-symmetrical 
static I — V dependence, measured at T = 1.285 K and H = 
—0.68 mT (close to the third matching field) for the sample 
A. The insets show the transformation of the sinusoidal bias 
current into a complex signal with nonzero average (the diode 
effect). 



and small microholcs ("antidots") of 1200 x 1200 nm^ 
and 600 x 600 nm^ in size (Fig. [T]). Depending on the 
separation s between centers of these antidots the result- 
ing complex array can be symmetrical (s = 1500 nm) 
or asymmetrical (s = 1100 nm) with respect to the 
y— axis. For brevity we introduced notations "sample 
S" and "sample A" (either symmetrical or asymmetrical 
pinning potential). 

Our ac four-probe measurements were carried out as 
follows: a sinusoidal driving current / = /g sin(27r/f), 
applied to such a bridge (the amplitude Iq ~ 141 [lK 
and the frequency / = 1.11 kHz), generates a voltage 
drop V{€). The typical 1 — V dependence is shown in 
Fig. [21 This time-dependent voltage can be represented 

oc 

as Vify = X) ^/ sin(27rri/t+iy9„), where Vnf and ipn are 

the amplitude and the phase of the n— th Fourier com- 
ponent. The amplitudes of the first, second and third 
Fourier harmonics V^/ , V2/ and Vij/ as well as the dc volt- 
age Vb as functions of the transverse magnetic field iJ 
and temperature T were measured by the lock-in ampli- 
fier Signal Recovery 7225 and the nanovoltmeter Hewlett 
Packard 34420A. All experimental results are summa- 
rized in Figs. [3] and [H 

We begin with the discussion of the H — T diagrams 
of the linear response and 1/^, measured for samples 
S and A at the fundamental frequency [panels (bl) and 
(b2) in Fig. [3] . The occurrence of resistivity minima at 
integer occupation numbers m = H/Hi in both samples 
indicates the formation of very stable vortex configura- 
tions. Here i/i = $o/a^ — 0.22 mT is the first matching 
field, at which the sample is threaded by one magnetic 
fiux quantum, <i>o = 2.07 • 10"^^ Wb, per unit ceU. The 



broader resistive transition of the sample S in compar- 
ison with the sample A (at the same T/Tc ratio) can 
be attributed to the enhanced surface superconductivity, 
since the antidots in the sample S are distributed more 
uniformly across the sample area. 

As expected, the loss of mirror symmetry of the pin- 
ning landscape along the axis significantly increases 
the dc rectification for sample A [compare the panels (al) 
and (a2) in Fig. [3]. The voltage pattern in the H ~ T 
diagram demonstrates the series of sign inversions which 
were recently explained in Ref. [lol | by taking into ac- 
count the intervortex interaction for a dense vortex lat- 
tice in the presence of a ratchet potential. Under the 
same conditions the diode effect for sample S is practi- 
cally absent, except in the vicinity of the third matching 
field {\H\/Hi c^i 3). The origin of such 'resonant' rectifi- 
cation for the nominally symmetrical sample is currently 
unclear. 

The H — T diagrams of the response at the second 
Fourier harmonic, 1^2/ ^-nd V^j seem to be similar to that 
for the dc response in many respects [compare the panels 
(al) and (a2) in Figs. [3] and [4]. Due to a larger "signal- 
to-noise" ratio for measurement of the second Fourier 
harmonic in comparison to dc measurements, this mod- 
ulation measurements appear to be more effective for a 
probing the symmetry properties of samples. Finally, we 
showed the results of the measurements at the triple fre- 
quency, Vif{H,T) and V^f{H,T) [panels (bl) and (b2) 
in Fig. H] . We emphasize that the field- induced oscilla- 
tions of the critical temperature, related to the formation 
of stable vortex structures, become more pronounced as 
compared with the Tc oscillations observed in the linear 
response. Thus, the experimental treatment of the non- 
linear response is a rather sensitive and powerful tech- 
nique for studying the commensurability effects in nanos- 
tructured superconductors. 

In order to illustrate correlations between the sym- 
metry of the pinning potential and the spectral com- 
position of the induced voltage, we simulated the low- 
frequency dynamics of a vortex chain in an asymmetrical 
ID periodic potential [lo| under the sinusoidal current 
excitation and estimated the instant mean vortex veloc- 

N 

ity {uy{t)) — X Vi/N {N is the number of vortices). 

i=l 

We considered both long-range vortex-vortex interac- 
tion, V{r) = —Eq ln(r/A), and short-range one, V{r) = 
Eq Ko{r/X) {Kq is the modified Bessel function), used re- 
cently for describing ratchet phenomena in superconduct- 
ing systems with the period comparable with the penetra- 
tion length A [l^ ■ The Fourier components of the instant 



mean vortex velocity. 



iy/nf 



"1///, 



as well as its dc value, (uy)o = /q {uy[t))dt, are ex- 
pected to determine the main contribution to the spec- 
trum of the voltage: Vb cx; H{uy)o and Vnf (x H(uy)nf- It 
was found that independently on the type of the intervor- 
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FIG. 3: (Color online) The dependencies of the voltage Vnf , measured at zero and fundamental on the external magnetic field 
H and temperature T: (al) Vq^ , (a2) Vq^, (bl) Vf , (b2) Vf. Solid black lines correspond to the conditions \Vnf{H,T)\ = 
a max V„f, a = 0.1, 0.5, 0.95. All insets show the dependencies of the Fourier components of the mean instant velocity (uy)„f 
and {Uy)nf on H and the effective inverse temperature U/Eq, calculated with the assumption of the logarithmic intervortex 
interaction V{r) = —Eo ln(r/A), where U is the pinning energy. Note that the periodic changes in {u^)o, (u^) f and {Uy) f at 
sweeping H resemble the variations of the Fourier components of the voltage. Color scales for all plots and insets are chosen 
in the same way. 



tex interaction, this simple model does explain how an ac 
excitation of interacting vortices in the ratchet potential 
leads to a generation of even Fourier harmonics observed 
experimentally. The results of our simulation are plot- 
ted in the H/Hi — U/Eq plane (see the insets in Figs. [3] 
HI, where U is the characteristic pinning strength. Since 
Up is known to decrease at increasing temperature and 
vanish at T = Tc we can consider a wide range of 
the pinning strengths (and temperatures), starting from 
Up = 0. It should be noted that our ID model cannot 
describe the rectification at \H\ ~ 37?i for the sample S. 
This indicates out that the dynamics of 2D vortex lattice 
in the 2D pinning potential is crucial for the explana- 
tion of all observed peculiarities in nonlinear response. 
However this issue goes beyond the subject of this Let- 
ter, and the detailed treatment in the framework of a full 



2D model similar to that used in Refs. [y, ll5| will be 
considered separately. 

Summarizing, we presented a comparative study of 
the linear and nonlinear properties of the perforated su- 
perconducting films. We demonstrated that the non- 
linear properties of the nanostructured superconductors 
can be tuned by making the pinning landscape of a re- 
quired symmetry, and thus we proposed a tunable vor- 
tex frequency converter. Interestingly, such a symmetry- 
breaking induced harmonics generation could be a basis 
for a designing the materials with an anisotropic nonlin- 
earity. 
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Fund GOA/2004/02 program, NFS - FSF program, the 
Belgian lAP, the Fund for Scientific Research - Flan- 
ders (F.W.O.-Vlaandcrcn), the Brazilian research agency 
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FIG. 4: (Color online) The dependencies of the voltage Vnf, measured at double and triple frequencies, on the external magnetic 
field H and temperature T: (al) Vj/, (a2) V^f, (bl) Vg^, (b2) Vg-). All notations are the same as in Fig. [3] 
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